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Abstract 

We show that we are able to extract the pion electromagnetic form factor by using the 
recent charged pion electroproduction data from JLab and a simple i-channel diagram. For 
this purpose we have performed the Q 2 -independent and dependent analyses. The result of 
the first analysis is in good agreement with those of previous works and fit the Maris and 
Tandy model as well as the monopole parameterization which describes a pion radius of 0.672 
fm. The result of the second analysis corroborates the findings in the first analysis. Our 
findings therefore provide a direct proof that at the given kinematics the i-channel diagram 
really dominates the process. This could also set a new constraint to the phenomenological 
models that try to describe the process. 
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1 Introduction 



In nuclear physics pions occupy a special place. Being the lightest meson in its family, pions play 
the major role in strong interaction. In fact, pions were the first mesons predicted by Yukawa as 
the mediator of the strong nuclear force. However, despite intensive investigations for decades, 
our understanding of the pion structure is far from complete. This is reflected, e.g., by the 
incomplete available data base of the pion electromagnetic form factor F n (Q 2 ) in a considerably 
wide range of virtual photon momentum transfer Q 2 . We understand that F w is nothing but a 
Fourier transformation of the charge distribution. 

At very small Q 2 the behavior of F % has been accurately determined by scattering of high 
energy pions from atomic electrons [Tj. At higher Q 2 values one should consult pion electropro- 
duction off a nucleon. 

With the continuous and excellent electron beams available at modern accelerators, such 
as CEBAF at JLab, it is now possible to measure pion electroproduction observables with 
unprecedented accuracy, at certain kinematics inaccessible with previous technology. This allows 
for a precise extraction of the longitudinal differential cross section via a Rosenbluth separation, 
at very small values of Mandelstam variable —t [21 [SJ 0j. These recent data show that the 
cross section do/dt at this kinematics exhibit a steep decreasing function of — t, indicating the 
dominance of the t-channel contribution in the process. This behavior becomes the first pillar for 
the F n extraction in Refs. [21 0. S] - The extracted F n values within the range of 0.6 < Q 2 < 2.5 
GeV 2 are impressively accurate. Nevertheless, it is important to note that these extractions: 

• rely heavily on the t-channel dominant assumption, although the assumption has never 
been directly checked. 

• are model dependent. In Refs. [21 El E] the Regge model of Ref. [5] has been utilized to 
obtain F n . 

• are plagued with the model uncertainties, which come from the extrapolation of the form 
factor cut-off to the physical limit. These uncertainties have been reported by Ref. [4j. 

• are "strategy" dependent. As shown in Ref. [4j, by using the same data and model as 
in Ref. [2] but focusing the analysis on the minimal t allowed by kinematics, where the 
unwanted background is assumed to be negligible, smaller F n values are obtained. 

The above facts indicate that current extraction methods are far from established. Since a 
completely model-independent extraction is quite difficult and challenging (although it is not 
impossible, see Ref. [6] for instance), attempts in this direction would be invaluable for helping 
to resolve this issue. 

It is the purpose of this paper to revisit a simple method, a la Chew-Low, which is based 
only on the first assumption. This will also serve as an alternative method, which is obviously 
easy to check. 

This paper is organized as follows. In Section [2] we briefly describe the formalism used in 
our calculation. Section [3] presents the result of the Q 2 -independent analysis, where we extract 
the electromagnetic form factor for fixed Q 2 values. Section 0J demonstrates the result of the 
Q 2 -dependent analysis, i.e., the result when we parameterize the form factor with a monopole 
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Figure 1: (Color online) The i-channel diagram that dominates the contributions to the pion electropro- 
duction process. The electromagnetic and hadronic form factors are denoted by and i*h, respectively. 



form and fit all available data simultaneously. We shall summarize and conclude our findings in 
Section [5l 



2 Formalism 

We begin with the convention of four-momenta given in the process 

lv (k)+N(p)^7T(q)+N(p') , 
to write the t-channel amplitude for the diagram shown in Fig. [T] as [TJ [H] 



(1) 



M 



2ieg nNN F n (k 2 )F h (t) 
k 2 (t — m 2 ) 



k ■ qk ■ e)u(p), 



where e is the proton charge, e denotes the virtual photon polarization, while m n and g % NN are 
the pion mass and the pion-nucleon coupling constant, respectively. 

In Eq. ([2]) the hadronic and electromagnetic form factors are denoted by F^ and F n , respec- 
tively, while the Mandelstam variable t can be written as 



t = (k-qf 



k 



2k q . 



(3) 



Note that in Eq. ([2]) we have added the so-called the Fubini-Nambu-Wataghin term [8j in order 
to retain the gauge invariance [9]. 

Clearly, this amplitude contributes only to the longitudinal cross section. For meson elec- 
troproduction it is well known that this cross section can be written as [10] 



da l -2\q\Wk 2 



(s — m 2 )/cQ 



l/sl 2 + l/el 2 + 2 cos ^Re(/ 5 */ 6 



where W indicates the total cm. energy, 9 is the pion scattering angle, and 
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with m denotes the nucleon mass and A§ is given by the term inside the square bracket in the 
r.h.s. of Eq. ([2]). Inserting Eq. ([5]) into Eq. ^ and by making use of 



(6) 



cLcfl vr d<JL 

~~dT ~ \%\ \q\ ~dM ' 

we obtain 

dt 87rW(s-m 2 )(t-m 2 ) 2 \k\ 3 V k 2 J 

Since we are working with small t, we can neglect both t and m n in Eq. © to obtain k 2 ~ 2k- q, 
so that the last term of Eq. (|7|) can be reduced to (go — |^o)- This has an advantage, because we 
have reduced the dependencies of the cross section on t. Finally, we can shift the pion propagator 
term to the l.h.s. of Eq. ([7]) and using Q 2 = —k 2 to obtain 

dt S-nW(s-m 2 )\k\ 3 K 2 ' 

The r.h.s. of Eq. (JSJ) were a linear function of t, if we could neglect the hadronic form factor 
Ffr. If this were the case, then an extrapolation of Eq. (JSj) to the pion pole (t = m 2 ) could be 
directly performed to determine the F w at a certain Q 2 value, because the r.h.s. is free of poles 
and the extrapolation (setting t at m 2 ) will eliminate other contributions beyond the t-channel 
one. However, as we clearly understand, this is not the case. 

To overcome the shortcoming of the extrapolation method, here we only assume that Eq. ([8|) 
is valid in the physical region (t < 0, but still small) because the t-channel dominates the process. 
Thus, for each value of Q 2 , Eq. © can be fitted to the data in the form of 

y = F*-f(t)+g, (9) 

to obtain F n . In Eq. Q we introduce the g parameter to account for other terms neglected 
in Eq. ([8]). Obviously, the use of this parameter improves the x 2 significantly, thus, absorbing 
more information from the electroproduction data, but at the cost of increasing the flexibility 
of the fit, i.e., increasing the error bars of the extracted F n . 



3 Q 2 -Independent Analysis 

In this analysis we first fix the value of the pion-nucleon coupling constant to be g^NN I '(^ 7r ) = 
14.17 [llj. We note that in the literature this value varies from 13.45 [TJ] to 14.52 [T3]- Thus, 
our choice lies nicely in this range. Later, we will also try to investigate the influence of the 
coupling constant variation on the extracted form factor. For the hadronic vertex we adopt the 
monopole form factor 

w) = { 1 -i 2 ) 1 • (10) 

Reference [T3| has pointed out that Ah = 0.85 GeV is generally used, whereas the precise value of 
Ah is a matter of controversy and may vary between 0.4 and 1.5 GeV. Meson exchange models for 
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Figure 2: (Color online) The modified longitudinal differential cross section as a function of the 
Mandelstam variable —t. Solid (dotted) lines show the results of fitting to experimental data with 
Ah = 0.85 (1.30) GeV. Solid squares [4] an d solid circles [3] show the recent JLab data, whereas solid 
triangles display the older data from DESY 19J. The upper (lower) numbers at the left part of each 
panel show the values of % 2 /N obtained for Ah = 0.85 (1.30) GeV. 
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Figure 3: (Color online) The electromagnetic form factors of pion extracted in this work and previous 
calculations. The monopole form factor is obtained by using the pion charge radius of 0.672 fm (PDG 
value 21). The dashed line shows the model of Maris and Tandy [22] . Results of the present work are 
obtained by using Ah = 0.85 GeV and shown in this figure by the solid (red) squares. Other data shown 
here are from Amendolia et al. pQ, Ackcrmann et al. [23j . Horn et al. [3], and Tadevosyan et al. 4J. The 
data point of Horn et al. at 1.6 GeV, as well as the corresponding result of the present analysis, have 
been slightly shifted for visual clarity. 



nucleon-nucleon interaction dictates that Ah > 1.30 GeV [15], whereas lattice QCD calculations 
prefer a value of 0.75 GeV [16) . On the other hand QCD sum rule calculations yield Ah = 0.80 
GeV [IT], while K. Vansyoc [18] found that Ah < 1.00 GeV is demanded by the present pion 
electroproduction data. 

In this work we found that the best result is obtained by using Ah = 1.30 GeV, i.e., the case 
when we consider the result of meson exchange models for the nucleon-nucleon interactions [15] . 
Nevertheless, using the generally accepted value (i.e. Ah = 0.85 GeV) the obtained result is still 
in good agreement with previous analyses. 

At W = 1.95 GeV there are four sets of data points from Ref. [I] with Q 2 = 0.6, 0.75, 
1.00 and 1.60 GeV 2 , as a function of t. At higher energy and momentum transfers (W = 2.22 
GeV, Q 2 = 1.6 and 2.45 GeV 2 ), two more data sets are available from Ref. [3]. For the sake of 
comparison, we also present the reanalyzed DESY data at W = 2.19 GeV, Q 2 = 0.70 GeV 2 [19] . 
These data were fitted separately by using the CERN-MINUIT code. A combination of the 
simplex and migrad packages has been used [20] . The fitted (modified) cross sections are shown 
in Fig. [21 where we can see that, except for the DESY data, the x 2 /N values obtained from the 
fits are quite encouraging. 

The extracted F w for Ah = 0.85 GeV and 1.30 GeV compared with the results of previous 
calculations are listed in Table Q] and depicted in Fig. [3j The error bars reported here are 
also obtained from MINUIT. Their magnitudes can be easily understood from the comparison 
between the fit results (solid and dashed lines) and the experimental data as demonstrated in 
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Table 1: Values of the extracted in the Q 2 -independent analysis obtained by using Ah = 0.85 GeV 
and 1.30 GeV compared with those of previous works. Reference [4] has reported "model uncertainties" 
listed as the second error-bars in the fourth column. 



Q 2 


Fir (this 


work) 


F w (previous works) 


(GeV 2 ) 


A h = 0.85 GeV 


A h = 1.30 GeV 




0.60 
0.70 
0.75 
1.00 
1.60 
1.60 
2.45 


0.470 ± 0.037 
0.467 ± 0.041 
0.435 ± 0.028 
0.393 ± 0.044 
0.251 ± 0.079 
0.337 ± 0.009 
0.209 ±0.017 


0.438 ± 0.066 
0.395 ± 0.076 
0.396 ± 0.088 
0.341 ± 0.074 
0.194 ±0.124 
0.276 ± 0.036 
0.153 ±0.042 


0.433 ± 0.017tolH 0! 
0.473 ± 0.023±g;g|| [4] 
0.341 ± 0.022 [H 
0.312 ± 0.016±g;gf| 0] 
0.233 ± 0.014lg;gjg 0] 
0.243 ±0.012 [3] 
0.167 ±0.010 [3] 



Fig. El 

Figure [3] reveals that, within the error bars, our simple analysis results in a good agreement 
with previous analyses [HJ [4] . At first glance the error bars obtained in the present work would 
seem to be much larger than those of previous works. However, this is in general not true. For 
example, consider the case of Q 2 = 0.6 GeV 2 . By combining the error bar obtained in Ref. [I] 
with the reported error bar coming from the "model uncertainty" we can clearly see that the 
error bar of the present work is indeed smaller. For other cases the error bars of the present 
work are mostly comparable to those of previous analysis. Only at Q 2 = 1.60 GeV 2 the obtained 
error bar is much larger. The reason behind this behavior is obvious from Fig. [2j Comparing 
with the result of Ref. [3] at the same kinematics we can conclude that this error bar can be 
significantly reduced if more statistics are available at this point. 

The two data points at 1.6 and 2.45 GeV 2 from Ref. [3] require a special explanation. First, 
Ref. [3] does not report the "model uncertainties" which appear during the extraction process. 
From Table Q] we can see that these uncertainties are substantial in increasing the error bars [4] . 
Second, the discrepancy appears between our result and those of Ref. [3] for these two data 
points are presumably due to higher W and t, where other terms excluded in Eq. ([8]) could 
start to contribute. Nevertheless, the error bars for the two F w at Q 2 = 1.6 GeV 2 of our results 
indicate that our analysis in this case is still consistent. Therefore, we believe that the error-bars 
obtained in the present work are realistic, especially if we compare them with those obtained 
from scattering of high energy pions from atomic electrons at Q 2 ~ 0.25 GeV 2 PQ. Thus, more 
statistics in the longitudinal cross section data will reduce these values. 

In Fig. 2] we demonstrate the variation of the extracted electromagnetic form factor F n due 
to the variation of the hadronic cut-off (Ah = 0.85 GeV and 1.30 GeV). The corresponding 
numerical values for the two cases are listed in Table [TJ A soft cut-off will substantially suppress 
the cross section. To overcome this suppression the fit will increase the electromagnetic form 
factor cut-off. This explains why the extracted F n values increase once we decrease the hadronic 
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Figure 4: (Color online) The influence of the variation of the hadronic form factor cut-off Ah on the 
extracted pion electromagnetic form factor. The monopole form factor is as in Fig. [3] 



cut-off. From this figure we may conclude that the best agreement with the results of Refs. [3] HJ 
is obtained if we set the hadronic cut-off to 1.3 GeV, a case when we consider the result of meson 
exchange models for the nucleon-nucleon interactions [15] . 

Note that in Figs. [3] and 0] we compare our results with the monopole form factor obtained 
by using the Particle Data Group (PDG) value [21] of the pion charge radius, i.e. 0.672 fm. In 
the literature one can find that the value ranges from 0.560 fm |24j to 0.740 fm [25] , whereas 
the recent Chiral Perturbation Theory extraction yields a value of 0.661 fm [26J. 

4 Q 2 -Dependent Analysis 

In the previous section we have extracted the electromagnetic pion form factors at fixed Q 2 
values. It is however possible to analyze the whole experimental data simultaneously in a fit. 
For this purpose, we have to assume that the electromagnetic form factor can be parametrized 
with a monopole form factor, 



A quick glance to Fig. 2] indicates that with this monopole parameterization it would be difficult 
to produce the pion charge radius as reported by PDG, especially if we used Ah = 0.85 GeV, for 
which we see that the two extracted form factors at Q 2 = 1.6 GeV 2 and 2.45 GeV 2 obviously 
overestimate the monopole parameterization. For Ah = 1.30 GeV we could expect that such a 
parameterization would still fairly work, since as shown in Fig. 2] the PDG form factor is still 
within the corresponding error bars. 

Nevertheless, a simultaneous fit to all available experimental data is quite important if we 
want to investigate the internal consistency of the data or the general trend of the extracted 
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Table 2: Values of some important parameters extracted from fit to three different sets of experimental 
data in the Q 2 -dependent analysis. References for the literature are given in the text. 



Parameter All data (36 points) 19 Data points 18 Data points Literature^ 

(r 2 ) 1 / 2 (fm) 0.576 ±0.020 0.596 ±0.111 0.625 ± 0.272 0.560 - 0.740 

g 0.019 ±0.006 0.022 ± 0.013 0.025 ± 0.007 

A h (GeV) 0.850 ± 0.034 0.968 ± 0.256 1.300 ± 0.230 0.400 - 1.500 

gl NN /(A-K) 14.500 ± 0.827 14.500 ± 0.623 14.500 ± 0.721 13.45 - 14.52 

X 2 /N 1.118 0.348 0.263 

^See text. 



parameters. Furthermore, within this framework we can leave the fit to determine not only the 
electromagnetic form factors F v , but also the hadronic cut-off Ah and the pion-nucleon coupling 
constant g n NN- Here, it is important to note that compared to the conventional approaches, 
which make use of thousands data points, a simultaneous fit to only 36 data points is clearly 
less reliable to accurately determine the hadronic cut-off and coupling constant of the pion. 
Therefore, our motivation here is only to confirm the results obtained by the Q 2 -independent 
analysis given in the previous section as well as to check our claim about the limitation of our 
model. 

As a first step, we fit all experimental data by leaving the pion charge radius (r 2 ) 1 / 2 (for 
the sake of comparison, instead of the A,,-), the constant g in Eq. ([9]), the hadronic cut-off Ah, 
and the coupling constant g-^NN, as free parameters. Since we only intent to qualitatively study 
the extracted parameters, we limit the pion charge radius in the range of 0.560 fm - 0.740 fm, 
the hadronic cut-off between 0.85 GeV and 1.30 GeV, and the pion-nucleon coupling constant 
9ttNN I between 13.50 and 14.50. The corresponding numerical values extracted from this 
fit along with their x 2, /N are shown in the second column of Table [2j The individual x 2 as a 
function of t is plotted in Fig. [5l Fitting to all available data results in a total \ 2 OI about 
36. Interestingly, we found that more than 36% of it come from the reanalyzed DESY data at 
W = 2.19 GeV, Q 2 = 0.70 GeV 2 [19]. This indicates that these data are probably incompatible 
with the recent JLab data [3j 0]. Furthermore, we found that these data belong also to the 
high-t and W data. 

The more surprising result is however shown by the set of the new JLab data at Q 2 = 2.45 
GeV 2 [3j. Although this data set contains five precise data points, the combination of high t, W, 
and Q 2 , obviously puts these data beyond the limitation of our approach. For Q 2 = 0.60 GeV 2 
only the data point at \t\ = 0.06 GeV 2 yields \ 2 > 1- This is understandable if we look at the 
first panel of Fig. [2J where it is obvious that this data point slightly underestimates the trend of 
the other data. The same situation happens for the last data point at Q 2 = 0.75 GeV 2 . Finally, 
we found that the first two data points at Q 2 = 1.00 GeV 2 and Q 2 = 1.60 GeV 2 overestimate 
the fit result. This fact is relatively difficult to explain, especially for the first case, for which 
the corresponding t are clearly small. Probably, these data points are too high or other physics 
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Figure 5: (Color online) Values of the % 2 obtained from fit to all experimental data in the Q 2_ dependent 
analysis plotted as a function of Note that the two numbers given in each panel show the corresponding 
values of Q 2 and W, respectively. The order of the experimental data sets shown in this figure (from top 
to bottom) is as in Fig. [2 



beyond our approach at this kinematics are missing. 

As a second step we remove the data points that produce > I from our fitting data 
base. The number of data points at this step turns out to be 19. The extracted parameters 
obtained from fitting to these data are shown in the third column of Table [2j Obviously the 
X 2 /N decreases significantly. The pion charge radius and hadronic cut-off increase slightly from 
those of the first fit, whereas the pion coupling stays the same. The only interesting result at 
this step is the fact that the highest x 2 comes from the data point at W = 2.22 GeV, Q 2 = 1.60, 
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Figure 6: (Color online) Electromagnetic form factors of the pion according to the results of fitting to 
three different sets of experimental data in the Q 2 -dependent analysis. The corresponding pion charge 
radii are listed in Table [2l For comparison, results of the Q 2 -independent analysis with Ah = 0.85 GeV 
and 1.30 GeV are also shown by the solid and open circles, respectively. 



\t\ = 0.22 GeV 2 . Again, this is understandable because of the combination of high t and W and, 
furthermore, this data point seems to underestimate the trend of other four data points. 

As the last step we remove the data point at W = 2.22 GeV, Q 2 = 1.60, \t\ = 0.22 GeV 2 from 
our fitting data base and refit the rest of the data (18 points). The extracted parameters for 
this step are shown in the fourth column of Tabled The value of X 1 /N is very small, indicating 
that our approach works pretty well for these selected data. The pion charge radius (as well as 
its error bar) and the hadronic cut-off increase significantly after we have removed the last data 
point as described above. Although having a relatively large error bar, the pion charge radius 
is now much closer to the PDG value compared to those obtained in the first and second steps. 
Thus, the result of the third step corroborates the finding in the Q 2 -independent analysis, i.e., 
smaller x 2 values would be obtained if we used A^ = 1.30 GeV. Another way to understand 
this result is by comparing the electromagnetic form factors F n obtained from this analysis with 
that of the PDG and those obtained from the (^-independent analysis described in the previous 
section. This is shown in Fig. [6j where we can clearly see that omitting the high-x 2 data (put 
in other words, using experimental data with the kinematics allowed by the model) results in a 
form factor that is in closer agreement with the form factor of PDG and those obtained from 
the Q 2 -independent analysis with Ah = 1.30 GeV. 

It is also important to note that the values of g and pion coupling constant are relatively 
stable for the three different fits (see Table [2]) . The relatively large error bars for the latter 
indicate that our conclusion drawn from the two analyses is not too sensitive to the variation of 
the g n NN within the fit limit. 
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5 Conclusion 



In conclusion, we would like to say that we are able to extract the electromagnetic form factor 
of pion from the recent longitudinal cross section data of pion electroproduction from JLab by 
using a simple i-channel diagram, which utilizes fewer assumptions compared to the previous 
calculations. For this purpose we have performed a Q 2 -independent analysis supplemented with 
a Q 2 -dependent analysis. The result of the first analysis shows a good agreement with F n 
extracted from previous works. The extracted values depend on the choice of the hadronic form 
factor cut-off. Nevertheless, by using the generally accepted value we still obtain a satisfactory 
result. The best result would be obtained if we used a cut-off from the meson exchange model. 
The result of the second analysis corroborates the findings of the first one. More accurate, with 
also smaller \t\, longitudinal cross section data will certainly help to reduce the uncertainty of 
the present calculation. Finally, we would like to mention that this simple calculation provides 
for the first time a direct proof that at the given kinematics pion electroproduction is dominated 
by the t-channel process. Previous works considered this merely as an assumption. This result 
could set a new constraint on the phenomenological models that try to explain the process, i.e. 
at the kinematics given in this paper contributions from other channels should vanish or a least 
should be minimal. We note that, however, at small t but higher energies, contributions of 
certain resonances to the process could be substantial. 
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